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Abstract 

It is observed that interstitial hydrogen nucleii on a metallic lattice are strongly 
coupled to their near neighbours by the unscreened electromagnetic field mediating 
transitions between low-lying states. It is shown that the dominant interaction is of 
dipole-dipole character. By means of numerical calculations based upon published 
data, it is then shown that in stoichiometric PdD, in which essentially all interstitial 
sites are occupied by a deuteron, certain specific superpositions of many-site product 
states exist that are lower in energy than the single-site ground state, suggesting the 
existence of a new low temperature phase. Finally, the modified behaviour of the 
two-particle wavefunction at small separations is investigated and prelimary results 
suggesting a radical narrowing of the effective Coulomb barrier are presented. 
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1 Introduction 



Several metallic elements, notably palladium, vanadium, niobium and nickel, 
can reversibly absorb hydrogen up to the point of stoichiometry, in which every 
available interstitial site - of octahedral (O) or tetrahedral (T) symmetry - is 
occupied by a hydrogen nucleus. A single hydrogen nucleus in such an environ- 
ment exhibits a spectrum of singlet, doublet and triplet state representations 
of the local point symmetry group. The ground state is invariably a singlet 
with even [+ + +] parity along each of the symmetry axes. The next level is, 

in an fee lattice, a triplet of states with parities [ — h +], [H h], [+ H — ] 

and an excitation energy of the order of 60 meV. The dipole moment between 
the ground state singlet and the first excited triplet is typically of the order 
of 0.2A Since the electronic Fermi gas couples weakly to the electromagnetic 
field quanta in this part of the spectrum, such a dipole moment gives rise to 
an essentially unscreened resonant dipole-dipole interaction (RDDI) between 



Preprint submitted to Elsevier 



4 February 2008 



nearest neighbours [lj. Simple geometrical considerations reveal this to be of 
the order of 20 meV per pair. Not only is this typically several times larger 
than screened static Coulomb interaction, it is also manifestly an appreciable 
fraction of the (on-diagonal) excitation energy of the dipole itself. Since in the 
quasi-stoichiometric loading regime each hydrogen has several nearest neigh- 
bours, it has previously been speculated [2] that there exist many-site states 
for which the total collective effect of the interaction is a multiple of the pair 
interaction. This paper sets out to answer the, in our opinion, intriguing ques- 
tion as to whether there exists any such collective state of quantum-entangled 
dipoles whose total energy is lower than the simple product of ground states, 
and to obtain an upper bound on the lowest possible energy of such an en- 
semble of coupled oscillators. 



2 Model 

The single particle states ip n are the solutions of 

ip n (r) = e n i> n (r) (1) 



H(r)4> n (r) 



2m 



V 2 + V(r) 



-where V(r) is the effective periodic lattice potential seen by a hydrogen nu- 
cleus. In view of the relatively large mass of hydrogen nucleii, the lowest energy 
solutions will generally be well-localised about local minima in V. These min- 
ima will coincide with sites of octahedral or tetrahedral symmetry in cubic 
lattices and hence the levels e n are an assortment of singlets, doublets and 
triplet representations of the cubic point symmetry groups. The static (zero 
frequency) components of the potential disturbance due to the hydrogen nu- 
cleus is subject to a screening law of the approximate (Thomas-Fermi) form 



-Kr 



V H h(t) 



(2) 



where K is proportional to the DOS at the Fermi surface. 
K is typically much greater than a reciprocal lattice vector. The static Coulomb 
H-H interaction between nearest neighbours is consequently small, typically 
not more than a few meV, and essentially state- independent. By contrast, the 
attenuation of the electromagnetic field due to a transition between levels is 
negligible over the dimensions of a lattice cell. 

Since the interparticle interaction is so strongly frequency dependent, the full 
Hamiltonian cannot be written in closed analytical form. However, matrix 
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elements between pairs of two-site states are simply given by: 

e -K\r 1 -r 2 -'R\S(e il -e^ )S(e i2 -ej 2 ) 

H h,h;i2,h = e 2 0'i , J 2 1 r TV (3) 

| r i — r 2 — rt| 

where R is the intersite displacement vector. The £(e» — e/) factors express the 
fact that only transitions between different levels give rise to an unscreened 
interaction. 

If just the lowest order term in the multipole expansion of the Coulomb op- 
erator 1 — 5-r is retained, there is no need for double integration over both 

|ri-r 2 -R| ' to 

volumes. In this approximation, ([3]) reduces to the familiar expression for a 
dipole-dipole interaction: 

e 2 3e 2 
H h,h;i2,j 2 ~ "BaO'ikiNi) • (J2|r 2 |i2> - ^(ji|R ■ ri|ii)0' 2 |R • r 2 |i 2 > (4) 



If a classical dipole is located at every interstitial O-site in an fee lattice, the 
interaction energy is lowest with the following orientations over a constant-z 
plaquette: 

°t !°t ! 

— > o < — o — > o < — o — > (5) 
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- where the open circles represent the sites of the metal cores at locations 
[100], [300], [111] etc. 

It can be shown that there is zero net interaction between parallel layers when 
each layer has such an arrangement. Guided by this classical analogue, we will 
limit our search for minimum energy states to those constructed from: 

[+ + +] parity states at all O-sites in a z = plaquette 
[ — h +] parity states at O-sites of even y and 
[H h] states at O-sites of odd y. 

[+ + +] ground states at all sites external to the plaquette 

For the rest of this paper we will use the shorthand \s,n) to denote the nth 
state of [+ + +] parity, \p x , n) to denote the rath state of [ — h +] parity and 

\p y ,n) to denote the nth state of [H h] parity. The singlet ground state is 

accordingly written as \s, 0). If the site location needs to be made explicit, we 
will append this in bold type thus: \s, 0, 110). For clarification, we reproduce 
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below an example of a pair of five-O-site states that are degenerate in zeroeth 
order and that are linked by the dipole-dipole interaction of (J4J: 



\s,0) 



\s,0) 



\Pv,0) 



\s,0) o \s,0) 
or equivalently: 



|s,0) o \ s ,0) 

o | s, 0) o 
\p*,0) ° |s,0) 



|s, 0, 000) g> |s, 0, 020) <g> |p„, 0, 110) ® |s, 0, 200) g> |a, 0, 220), 
1^, 0, 000) ® |s, 0, 020) ® |s, 0, 110) ® |s, 0, 200) ® [s, 0, 220) 

The Hamiltonian matrix in the subspace of the two-O-site states 

\s, 0,000) ® |p y , 0,110) and 

jpai, 0, 000) ® |s, 0, 110) is, according to (gj): 



/ 



H 



3V2e 2 d 2 



e s,0 + e p,0 
^3 e s,0 -T 6p,0 



(6) 



(7) 



where the dipole moment d = (s,0\x\p x ,0) = (s,0\y\p y ,0) 

The energy eigenvalues are in this case simply e Sj o + e Pi o ± 3y ^ c f d2 

More generally, and assuming full hydrogen occupancy and perfect lattice 

symmetry, the full many-site Hamiltonian matrix depends numerically upon 

just the lattice parameter a, the dipole moments d and the energies e of the 

single-site states. 



3 Application to PdD 



Quasi-stoichiometric PdH and PdD are natural candidates for our model be- 
cause the adiabatic effective potential experienced by the hydrogen nucleus 
has been determined from ab initio DFT calculations [3|4] and the theoretical 
spectra found to agree well with IR spectroscopic measurements over a wide 
range of substoichiometric loading ratios. For this work we checked the results 
published in [5] by solving (jTJ) on a real-space wedge mesh of pitch 0.03 A, us- 
ing the published effective adiabatic Pd-H potential and boundary conditions 
appropriate to the desired parity. The lattice parameter, corresponding to a 
displacement like [200] in our notation, is 4.07 A. The lowest few eigenvectors 
of the very large sparse matrix were solved using the dominant-diagonal iter- 
ation method, as in [5]. 
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The lowest single-site energy levels, relative to the O-site potential, were found 
to be: 



level 


PdH (meV) 


PdD (meV) 


e s,0 


82 


52 


e p,o 


151 


95 




233 


149 




289 


186 



- with the following dipole moments along each of the three cartesian axes: 



dipole 


PdH (A) 


PdD (A) 


(* 


0\x\p x ,0) 


0.172 


0.153 


(s 


l\x\p x ,0) 


0.128 


0.112 


(s 


0\x\p x , 1) 


-0.003 


-0.002 




l\x\p x , 1) 


0.194 


0.165 



Double integration over both site volumes according to (j3J) yielded the follow- 
ing off- diagonal elements: 
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^»i,ji;»2,j2 


PdH (meV) 


PdD (meV) 


(s 


0,000 


<8> (Pr 





110\H\p v , 


0,000) 




5,0,110) 


-27 


-21 


(s 


0,000 


<g) (p x 





110\H\p v , 

\ \± y ' 


0,000) 


<g> 


s, 1,110) 


-20 


-15 


(s 


1,000 


<£> (Px 





110\H\Py, 


0,000) 


<s> 


s, 1,110) 


-15 


-11 


(s 


0,000 


<E> (Px 


1 


U0\H\ Py , 


0,000) 


<s> 


s, 1,110) 


-30 


-23 


(s 


0,000 


<8> (Pr 





2001^1^ 


0, 000) 


<s> 


s, 0,200) 


-13 


-10 


(s 


0,000 


<8> (Px 





200|#|p x 

1 1 1 •>- 


0, 000) 




s, 1,200) 


-9 


-7 


(s 


1,000 


® (Px 





200\H\p x 

1 u J- 


0, 000) 




s,l,200) 


-7 


-5 


(s 


0,000 


® (Px 


1 


200|#|p x 


0, 000) 




s, 1,200) 


-14 


-11 


(s 


0,000 


® (p x 





020|#|p x 


0, 000) 




s, 0,020) 


6 


5 


(s 


0,000 


® (p x 





020|F|p x 

1 1 1 •>- 


0, 000) 




s,l,020) 


5 


4 


(s 


1,000 


® (j) x 





020|#|p x 

1 1 1 •>- 


0, 000) 




s, 1,020) 


4 


3 


(s 


0,000 


® (Pr 


1 


020\H\p x 


0,000) 




s, 1,020) 


7 


5 


(s 


0,000 


® (Px 





220|#|p x 


0, 000) 


® 


s, 0,220) 


-1 


-1 


(s 


0,000 


<g> (p x 





220|#|p x 


0, 000) 




s,l,220) 


-1 


-1 


(s 


1,000 


<S> (]?x 





220|#|p x 


0, 000) 




s,l,220) 


-1 


-1 


(s 


0,000 


® (Px 


1 


220|#|p x 


0, 000) 




5,1,220) 





-1 



3.1 5-site states 



As an illustrative example, we will calculate the lowest 5-site energy achievable 
for a given number of \p, 0) states. There are just five 5-site states having 1 
\p, 0) and 4 |s,0) states, namely: 



\s, 0, 000) <g> \s, 0, 020) <g> |s, 0, 110) ® \s, 0, 200) <g> |p x , 0, 220) 
|s, 0, 000) <g> |s, 0, 020) (g) |s, 0, 110) <g> \p x , 0, 200) <g> \s, 0, 220) 
|s, 0, 000) <g> |s, 0, 020) <g> |p y , 0, 110) <g> |s, 0, 200) <g> |s, 0, 220) 
|s, 0, 000) <g> |p x , 0, 020) <g> \s, 0, 110) <g> |s, 0, 200) <g> |s, 0, 220) 
\p x , 0, 000) <g> \s, 0, 020) <g> |s, 0, 110) <g> |s, 0, 200) <g> |s, 0, 220) 



(8) 
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The corresponding Hamiltonian matrix is (in meV): 



43 5 -21 -10 -1 
5 43 21 -1 -10 
H= -21 21 43 21 -21 
-10 -1 21 43 5 
-1 -10 -21 5 43 



\ 



5e s 0I5 



(9) 



/ 



- where we have subtracted out the energy of the conventional ground state. 
The energy eigenvalues in this small subspace are 3,29,37,59 and 88 meV 
relative to 5e S)0 . For the 80 five-site states that comprise just one \p, 0) and 
four \s,n < 2) states, the lowest two energy eigenvalues are found to be -2 
and 26 meV relative to 5e S)0 . It is clear from this that the \s, 1) states make a 
significant contribution to the lowest energy eigenvector. 
In an attempt to find a converged value for the absolute mimimum site energy, 
Hamiltonian matrices were constructed for a series of plaquettes of increasing 
size up to a limit set by the memory capacity of our machine. The 15 sites 
included were: [000], [020], [110], [200], [220], [130], [310], [330], [420], 
[1-10], [-110], [240], [130], [040], [150] , in that order. 
The results obtained for the energies (Eq, E\) of the lowest two states relative 
to iVe S) o are summarized in the following table. For the larger plaquettes, an 
energy cut-off was applied in order to limit the size of the matrix. 
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Sites 


p-states 


Cut-off (meV) 


States 


E (meV) 


E 1 (meV) 


E / Site (meV) 


5 


1 




80 


-2 


26 





5 


2 




80 


31 


43 


6 


6 


1 




192 


-4 


19 


-1 


6 


2 




240 


12 


37 


2 


8 


1 




1024 


-15 


12 


-2 


8 


2 




1792 


-9 


12 


-1 


9 


1 




2304 


-16 


9 


-2 


9 


2 


300 


1044 


-14 


5 


-2 


12 


1 


300 


2784 


-21 


-7 


-2 


12 


2 


300 


3696 


-33 


-9 


-3 


12 


3 


300 


2200 


-34 


-12 


-3 


12 


4 


300 


4455 


-24 


-8 


-2 


15 


2 


300 


9660 


-40 


-23 


-3 


15 


3 


300 


5915 


-44 


-25 


-3 


15 


4 


250 


1365 


-26 


-8 


-2 



4 Conclusions 



The intrinsic complexity of this exact method and the inapplicablity of a per- 
turbative approach have so far confounded our attempts to establish a lower 
bound on the absolute minimum site energy. It follows from the variational 
principle that inclusion of higher \s, n) states, as well as further increase in pla- 
quette size, will result in even lower minimum energies. A mean-field approach 
is perhaps indicated, but we have as yet to find a sufficiently accurate formula- 
tion. It is nevertheless already clear from the above data that entangled states 
are favoured in the stoichiometric regime. The existence of a low temperature 
phase in which all the deuterons cohere in a mesoscopically entangled state is 
hence strongly indicated. 
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5 Further work 



At small interparticle distances :- K\rx — r 2 \ << 1, the off-diagonal elements of 
the type we have been considering are comparable in magnitude, but opposite 
in sign, to the static Coulomb pair repulsion term. It is hence reasonable so 
suppose that, once coherence has been established, the height and width of 
the effective Coulomb barrier between neighbouring s,p state pairs is reduced, 
with a concommitant increase in the - normally infinitesimally slow - D-D 
fusion rate. In order to investigate the magnitude of this effect, we solved 
the two-particle Hamiltonian for the two states discussed in connection with 
(J7|) above. Both the static (j2j) and dynamic ([3]) interactions were included. 
Memory constraints limited us to a grid resolution of 0.06A In view of (J6]) 
and (|7j), the solution was constrained to be of the form: 



ip(x 1 ,y 1 ,z 1 ,x 2 - f, 2/2 - f,^ 2 ) 

K -4>(x 2 - f,y 2 - f,3i,yi,zi,*i) 




where ip is odd in its 5th argument and even in all others. 
It was found that the lowest energy solution, with 

3V2e 2 d 2 

e ~ e s ,o + e P ,o ; (H 

a 6 



- exhibited an increased probability for close encounters of the two hydrogen 
nucleii right down to the limit of our resolution. At \ri — r 2 \ = 0.06A, the am- 
plitude was enhanced by about an order of magnitude over the simple product 
state that pertains when interaction is neglected. The region of overlap was 
concentrated about the T-site lattice potential minima that are equidistant 
between the two O-sites. Accessing the interesting region |r*i — 7*2] < 4 fm, 
demands a multigrid approach that is currently in progress. 
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